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We discuss recent results for the phase transition in finite-temperature QCD from numerical (Monte Carlo) 
simulations of the lattice-regularized theory. Emphasis is given to the case of two degenerate light-quark flavors. 
The order of the transition in this case, which could have cosmological implications, has not yet been established. 

PACS numbers: ll.lO.Wx, 12.38.Aw, 12.38.Gc, 68.35.Rh 

Keywords: QCD phase transition, lattice formulation, numerical simulations, universality, scaling 



o 
o 

(N 

Oh' 
D 



> 

m • 
o 

O . 
\o • 
o 



X 



I. INTRODUCTION 

Although we are still trying to understand how confine- 
ment arises in QCD, we know that under normal conditions 
the running coupling is large at low energies, causing quark 
and gluons to be confined. However, when exposed to ex- 
treme conditions such as very high temperatures or densities, 
quarks are forced to stay at very short distances from one an- 
other and there is a transition to a deconfined, quark-gluon- 
plasma phase. This transition is present both at high tempera- 
ture and at high density, implying a phase diagram where the 
hadronic phase exists only near the origin of the plane defined 
by the temperature and chemical-potential axes. Of course, 
one wants to know the location of the transition in order to 
study properties of the new phase of matter, but it is also inter- 
esting to determine the nature of this transition. In particular, 
it is important to establish if the transition is a strong one, of 
first order, involving discontinuity in the order parameter, or 
if it is such that the two phases are connected smoothly. This 
may have consequences for the understanding of the cosmo- 
logical QCD phase transition, which occuiTed a few microsec- 
onds after the big bang and formed the hadrons we observe 
today. In this case the transition lies closer to the tempera- 
ture axis and its nature is of direct importance to determine 
the types of cosmological rehcs that can be associated with it. 
In particular, a first-order transition would very likely be as- 
sociated with the formation of cold dark matter clumps JHI^l- 
The nature of the transition must be also taken into account 
at relativistic heavy-ion collision experiments. An additional 
requirement in this case is the description of dynamic effects 
il. 

The task of studying the QCD phase transition theoretically 
must be carried out by nonperturbative methods and a natural 
choice is to consider the lattice regularization as a formulation 
of the theory. In fact, lattice-QCD simulations allow a nonper- 
turbative description of the phase transition in hadronic matter 
at high temperatures and there has been some recent progress 
in the description of the transition also in the case of finite 
density 0|- In the case of the finite-temperature transition, 
there is a qualitative difference when dealing with the full- 
QCD case (i.e. considering dynamic quarks) or with the so- 
called quenched case, in which the gluonic effects are taken 
into account but sea quarks are taken to be infinitely massive. 



For the quenched case one studies the deconfining transition 
itself, by means of the order parameter given by the Polyakov 
loop. The transition in this case is of first order. In the fuU- 
QCD case there is no equivalent order parameter for the de- 
confinement transition and one must consider the chiral phase 
transition. This transition occurs when the chiral symmetry 
— exact in the limit of zero quark masses and spontaneously 
broken at low temperatures — is restored at high temperature. 
The case of two dynamic quarks, i.e. considering dynamic ef- 
fects of only two degenerate light-quark flavors, correspond- 
ing to the up and down quarks, is particularly interesting. In 
this case, if the transition is of second order, one would ex- 
pect to observe universal critical scaling in the class of the 
3d (9(4) continuous-spin model j5l|6|]. Also, in the contin- 
uum limit, simulations using different discretizations for the 
fermion fields should give the same results. The fact that the 
critical behavior should be in the universality class of a spin 
model can be precisely checked, since the nonperturbative be- 
havior for these models can be obtained with Monte Carlo 
simulations by so-called global methods, which avoid the crit- 
ical slowing-down present in QCD simulations 

The determination of the correct nature of the transition in 
the two-flavor case is one of the present challenges of lattice 
QCD, as pointed out by Wilczek in |9|. This prediction has 
been investigated numerically by lattice simulations for over 
ten years, yet there is still no agreement about the order of the 
transition or about its scaling properties H fioll . More pre- 
cisely, the predicted (9(4) scaling has been observed in the 
Wilson-fermion case but not in the staggered-fermion 
case, believed to be the appropriate formulation for studies 
of the chiral region. In this case, extensive numerical stud- 
ies and scaling tests have been done by the Bielefeld |12], 
JLQCD 1 13] and MILC Q groups. It was found that the 
chiral-susceptibility peaks scale reasonably well with the pre- 
dicted exponents, but no agreement is seen in a comparison 
with the (9(4) scaling function. At the same time, some re- 
cent numerical studies with staggered fermions claim that the 
deconfining transition may be of first order 

(Bill- 

In fnl ITsil a simple method was introduced to obtain a 
uniquely defined normalization of the QCD data, allowing 
an unambiguous comparison to the (normalized) (9(4) scaling 
function. The analysis showed a surprisingly better agreement 
for the larger values of the quark masses. Let us note that in 
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previous scaling tests the comparison had been done up to a 
(non-universal) normalization of the data and a match to the 
scaling function was tried by fitting it to the data points with 
the smallest masses. One interpretation of this result is that 
data at smaller masses (closer to the physical values) suffer 
more strongly from systematic errors in the simulations. In 
fact, larger quark masses are much easier to simulate, allow- 
ing greater control over errors and more reliable results. Here 
we present a preliminary study at a rather large mass value 
(nici = 0.075 in lattice units), using staggered fermions and 
the MILC code. We consider the standard action and tempo- 
ral lattice extent A'x = 4, as in the studies mentioned above. 

II. SCALING TESTS 

The behavior of systems around a second-order phase tran- 
sition (or critical point) may show striking similarities for 
systems that would otherwise seem completely different. In 
fact, it is possible to divide systems into so-called universality 
classes, in such a way that each class will have, e.g., the same 
critical exponents around the transition. Typical exponents are 

M,,=o.,^o- ^ kP, (1) 

Xh=o.t^o ^ kr^, (2) 

M,=o,/,^o ^ /^'/^ (3) 

where M is the order parameter — e.g. the magnetization for 
a spin system — % is the corresponding susceptibility and 

t = {T-T,)/To, (4) 
h = H/Ho (5) 

are the reduced temperature and magnetic field, respectively. 

Thus, in principle, one may compare the critical exponents 
for different systems to check if they belong to the same uni- 
versality class. In practice, however, the critical exponents 
may vary little from one class to the other and in order to carry 
out the comparison one would need to have a very precise de- 
termination of the exponents, which is not yet feasible in the 
QCD case. 

A more general comparison is obtained through the scaling 
functions for both systems. This comparison allows a more 
conclusive test, and can be applied for cases where the crit- 
ical exponents cannot be established with great accuracy. In 
this case we may assume the exponents for a given class and 
compare the behavior of the whole critical region for one sys- 
tem to the known scaling curve for the proposed universality 
class. The scaling Ansatz is written for the free energy Fs in 
the critical region as 

Fsit,h) = b-''F,iby'ty"h), (6) 

where b is a rescaling factor, d is the dimension and yt,yi, 
are related to the usual critical exponents: P, y, 5 mentioned 
above. The scaling Ansatz implies that the order parameter 
must be described by a universal function 

M//2'/5 = /^(f//,l/p5), 



The statement that the function /m is universal means that 
once the non-universal normalization constants 7b and Hq are 
determined for a given system in the universality class, the 
order parameter M scales according to the scaling function /m 
for all systems in this class. As said above, the comparison of 
(normalized) scaling functions between two systems is a more 
general test of universality, especially in the case of the QCD 
phase transition. 

A further difficulty in studying the critical behavior at the 
QCD phase transition is the impossibility of considering the 
critical point directly, since that would correspond to having 
zero quark mass, or zero magnetic field H in the language of 
the spin models above. In order to check scaling with critical 
exponents of a given class, or to determine the normalization 
constants To and Ho for systems where a study at // = is not 
possible, it is important to determine the pseudo-critical line, 
defined by the points where the susceptibility % shows a (fi- 
nite) peak. This corresponds to the rounding of the divergence 
that would be observed for H — Q, T — Tc. The susceptibility 
scales as 

X = dM/dH = {\/Ho)h'l'>-'f^{t/h"^) , (8) 

where /j^ is a universal function related to /m- At each fixed h 
the peak in % is given by 

tp = Zph'/^\ (9) 

Mp = h'^^fMizp), (10) 

HoXp = h'/^-'Mzp). (11) 

Thus, the behavior along the pseudo-critical line is determined 
by the universal constants Zp, fmizp), fxizp)- Critical ex- 
ponents, the scaling function /m and the universal constants 
above are well-known for the 3d (9(4) model 121.. 1221. .23.1 . 

Note that one may also consider the comparison for finite- 
size-scaling functions, since they are also universal and have 
the advantage of being valid for finite values of L, the linear 
size of the system. Such functions were studied for the 3d 
0(4) model in iQ. 

III. COMPARISON OF QCD DATA WITH THE 
PREDICTED SCALING FUNCTION 

We now turn to the comparison of the two-flavor QCD data 
in the critical region (in the case of small but nonzero quark 
mass) to the predicted scaling properties of the 3d 0(4) spin 
model. As mentioned in the Introduction, we consider the chi- 
ral phase transition, since there is no clear order parameter for 
the deconfinement transition in the case of full QCD. The or- 
der parameter for the chiral transition is given by the so-called 
chiral condensate < Tj7\(/ >, where i|/ is a combination of the 
quark fields entering the QCD Lagrangian |6]. The analogue 
of the magnetic field is the quark mass m^, and (on the lattice) 
the reduced temperature is proportional to 

6// -6/^2(0), (12) 

where g is the lattice bare coupling and g^ (0) is its extrapo- 
lated critical value. Therefore, referring to the pseudo-critical 
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line described in the previous section, the chiral susceptibihty 
peaks at 



h-^m,'l^\ (13) 



As mentioned in the Introduction, previous results from 
lattice-QCD simulations in the two-flavor case show good 
scaling (with the predicted exponents) only along the pseudo- 
critical line, which is given by the peaks of the chiral suscep- 
tibility. It should be clear from the discussion in the above 
sections that this is not a sufficient test to prove that the tran- 
sition is second order, especially if no agreement is seen when 
comparing the data to the scaling function. As described in 
we use the observed scaling along the pseudo-critical 
line and the universal quantities Zp, fM{zp) from the (9(4) 
model to determine the normalization constants Ho, Tq for the 
QCD data. This allows an unambiguous comparison of the 
data to the scaling function /m. More precisely, we note that in 
previous analyses 1 14] the normalization constants were tenta- 
tively adjusted by shifting the (9(4) curve so as to get a rough 
agreement with the data at smaller quark masses, since these 
are closer to the chiral limit. The problem is that the lighter 
masses are also more subject to the presence of systematic er- 
rors in the simulations. In this case the overall agreement was 
rather poor, indicating that there were strong systematic effets 
or that the transition is not in the predicted universality class. 
Here we fix the constants as described in Section |II| follow- 
ing the behavior along the pseudo-critical line. In this way no 
value of the quark mass is priviledged and the comparison is 
unambiguous. 

Our comparison is shown in Fig. below. The pseudo- 
critical line corresponds to a point in this plot and is marked 
with an arrow. For clarity we do not show the data (from the 
Bielefeld and JLQCD collaborations) obtained directly at the 
pseudo-critical point. These are slightly scattered around Zp 
but show good scaling within errors. We see relatively good 
scaling in the pseudo-critical region, i.e. around [zp, fi^izp)}, 
as expected. Away from this region most MILC points are 
several standard deviations away from the predicted curve. 
These data are given for three values of the quark mass in 
lattice units: 0.008, 0.0125 and 0.025. Note that the points 
with larger mass come closer to the curve. In particular, we 
can see that the new data at — 0.075 show sensibly bet- 
ter scaling, especially for larger temperatures, where previ- 
ously the scaling seemed unlikely. The good agreement of 
these data with the (9(4) scaling function motivates a care- 
ful study of systematic errors for smaller masses. A possible 
source of such errors are finite-size corrections, which would 
be stronger for smaller masses, since then the lattice side may 
not be large enough to "contain" the physical particle. Put dif- 
ferently, finite-size effects are expected when the correlation 
length (in lattice units) associated with a particle is compara- 
ble to or larger than the lattice side. Of course, this is more 
likely to occur for a lighter particle. 



IV. FINITE-SIZE EFFECTS 

In addition to the infinite-volume scaling laws mentioned 
above, we may also consider finite-size-scaling functions. In 
fact, the scaling Ansatz in Eq.|6]also implies 

M = L-P/^e,(/!LP^/^) (14) 

where L is the linear size of the system and we consider fixed 
values of the ratio t /h^l^^ = z (e.g. z = as in the critical 
isotherm, or Zp as along the pseudo-critical line). Thus, M can 
be described by a universal finite-size-scaling (FSS) function 
of one variable. We note that in order to recover the infinite- 
volume expression M = h^^^ fM{z) as L ^ we must have 
Qziu) fM{z)u^^^ for large u. Thus, in this limit, the FSS 
functions are given simply in terms of the scaling function 
/m(z)- Working with the FSS functions 2- instead of the 
infinite-volume scaling function /m has the disadvantage that 
one must consider z fixed (thus restricting the regions to be 
compared in parameter space) but the advantage that a com- 
parison can be made already at finite values of L. 

A finite-size-scaling analysis as described above was car- 
ried out in 1 17], but it was found that the QCD data show 
good (finite-size) scaling only along the pseudo-critical line. 

V. CONCLUSIONS 

Understanding the nature of the chiral phase transition in 
two-flavor QCD has proven to be a challenging task. The pre- 
diction of a second-order transition with critical behavior in 
the universality class of the (9(4) spin model is not verified 
for staggered fermions of small masses, although it can be 
shown (by an unambiguous normalization of the data) that 
better scaling is obtained for the existing data at larger (un- 
physical) masses. The fact that data for heavier quarks would 
show such good scaling was not expected, since the normal- 
ization of the data for comparison with the scaling curve did 
not priviledge any particular values of the quark mass. This 
suggests that the lack of scaling at small masses observed so 
far may be due to systematic effects, which could be due to 
finite-size corrections or to uncontrolled errors in the hybrid 
Monte Carlo algorithm used for updating the configurations. 
Both these sources of errors would be more significant for 
the case of smaller masses. As discussed above in Section 
IIVI above, the deviations observed are most likely not due to 
finite-size corrections and we thus argue that the deviations 
from (9(4) scaling at smaller masses may come from system- 
atic errors in the simulation, probably related to the use of the 
R algorithm for the simulations ]24]. Note that, contrary to 
what happens for the quenched-QCD case, the algorithm used 
to update full-QCD configurations is not exact and should 
have its accuracy tested carefully for each different value of 
the quark mass used. 

Let us also mention that a redefinition of the reduced tem- 
perature in terms of the physical temperature T including a 
term in the quark mass M, as suggested in ] 16], improves the 
agreement with the scaling curve further, as has been recently 
shown in L20il . 
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FIG. 1: Comparison of QCD (staggered) data to the 0(4) scaling function. For clarity, we do not show the data around the pseudo-critical 
point (indicated by the arrow), which were used to determine the normalization of the remaining data points. 
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